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ABSTRACT 

We report on the successful science verification phase of a new observing mode 
at the Keck interferometer, which provides a line-spread function width and sam- 
pling of 150 km/s at X'-band, at a current limiting magnitude of K' ~ 7 mag 
with spatial resolution of A / 2 B ps 2.7 mas and a measured differential phase sta- 
bility of unprecedented precision (3 mrad at K = 5 mag, which represents 3 //as 
on sky or a centroiding precision of 10 -3 ). The scientific potential of this mode 
is demonstrated by the presented observations of the circumstellar disk of the 
evolved Be-star 48 Lib. In addition to indirect methods such as multi-wavelength 
spectroscopy and polaritmetry, the here described spectro-interferometric as- 
trometry provides a new tool to directly constrain the radial density structure in 
the disk. We resolve for the first time several Pfund emission lines, in addition to 
Br 7 , in a single interferometric spectrum, and with adequate spatial and spec- 
tral resolution and precision to analyze the radial disk structure in 48 Lib. The 
data suggest that the continuum and P /-emission originates in significantly more 
compact regions, inside of the Br 7 emission zone. Thus, spectro-interferometric 
astrometry opens the opportunity to directly connect the different observed line 
profiles of Br 7 and Pfund in the total and correlated flux to different disk radii. 
The gravitational potential of a rotationally flattened Be star is expected to in- 
duce a one-armed density perturbation in the circumstellar disk. Such a slowly 
rotating disk oscillation has been used to explain the well known periodic V/R 
spectral profile variability in these stars, as well as the observed V/R cycle phase 
shifts between different disk emission lines. The differential line properties and 
linear constraints set by our data are consistent with theoretical models and lend 
direct support to the existence of a radius-dependent disk density perturbation. 
The data also shows decreasing gas rotation velocities at increasing stello-centric 
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radii as expected for Keplerian disk rotation, assumed by those models. 

Subject headings: stars: emission-line, Be — circumstellar matter — stars: individual 
(48 Lib) — techniques: interferometric — techniques: spectroscopic — infrared: stars 
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Introduction 



The ASTRA (ASTrometric and phase- Referencing Astronomy) upgrade program, 
funded by the NSF Major Research Instrumentation program, aims at extending the 
sensitivity and spectral resolution of the Keck interferometer (KI) through phase referencing 



and to im plement a narrow-angle 



of the KI ( Shao fc Colavita 



1992 



astrometry mode to enable a much 



Droader scientific use 



Wizinowich et al. 



2007 



Pott et al. 



2008 



Woillez et al. 



2010J). In this paper we report on results of the successful implementation of the first 
ASTRA mode, the self phase referencing (SPR) mode. In the SPR mode, 55 % of the light 
from the on-axis science target is used to provide fringe tracking (i.e., phase referencing) 
information while 20 % is used to make the science measurement. Due to the stabilized 
fringes, >100x longer integrations can be taken on the second camera. Therefore, it is 
possible in SPR to increase the spectral resolution of the interferometric measurement 
by at least an order of magnitude. SPR as currently implemented and provided to the 
general user, offers a spectral resolution R(= A/AA) of 2000, which is Nyquist sampled 
at R = 1000 (330 pixels across the K'-band) at K < 7 mag. Further tec hnical aspects 



Woillez et al. 



( 1201 Oh . Besides simple 



and results from the commissioning are discussed by 
amplitude spectro-interferometry at this resolution, the here presented data demonstrate 
that the differential visibility phase can be retrieved at a precision of at least 3 mrad for 
a bright star. Non-linear phase changes in the source spectrum can be measured, and 
typically indi cate translations o 



the photo-center over the respective spectral channels on 
20071 ). The operational readiness of KI-SPR was successfully 
demonstrated with the here discussed observations of 48 Lib on the night of April 25, 200 8 



the sky (e.g. 



Wei gelt et al. 



(UT) . In addition, several young stellar objects were observed successfully ( lEisner et al. 
20101 ) . We show at the example of 48 Lib the wealth of spatial and spectral information 
that can be provided by SPR data. 
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48 Lib is a well studied classical Be-star with circumstellar shell-emission lines. Previous 
optical-nearinfrared long-baseline interferometric (OLBI) measurements of Be-star-shells 
demonstrated that 3d-models for H-line and continuum emission often overpredict the 



true size, which demonstra tes the need for direct OLBI measurements (IStee fc Bittar 



2001 



Chesneau et al 



( Shao et al 



20051 ). Narrow-band observations with the Mark III interferometer 



1988), result ed in Be-star disk diameter estimates ranging in 2.6 to 4.5 mas 



( IQuirrenbach et al 



19971 ). well suited for the angular resolution of the 85 m baseline of the 



KI (X/2B = 2.7 mas). We chose 48 Lib as science demonstration target because recent 



detection of very narrow sat ellite absorption features in 



very high inclination angles ( jHanuschik fc Vrancken 



e II shell lines reliably indicate 



19961 ) . This ensures strong features in 



the differential visibility and phase signals across the emission lines. 

A large fraction of Be stars show cyclic variations in the ratio between the violet and 
red flux of the HI emiss ion lines (V/R), an enigmatic phenomenon studied since almost 



a century (ICurtiss 



19251 ). Early measurements with the GI2T interferometer spatially 



and timely resolved the Ha emission o 



emission region with time 



Mourard et al. 



7-Cas , showing spatial variation of the line 



(1989) 



rotational velocities for cooler stars like 48 Lib. 



. The fast rotation of Be-stars (near-critical 



Cranmer 



20051 ) r esults in a flattening of the 



stella r shape, which can be observationally confirmed by OLBI (jDomiciano de Souza et al. 



20031 ). While the rapid rotation might play a role in the ejection of the circumstellar 



envelope, the quadrupole moment of the gravitational potential of a flattened Be star is 



expected to create a prograde on e-armed 



spira. 



disk and creating the V/R cycles ( jOkazaki 



1991 



density pattern, 



slowly precessing in the 



Papaloizou et al 



1992). The latter authors 



calculate that the radial extension of such a disk oscillation mode should be confined to a 
few stellar radii, which can be probed by interferometry. Central properties of this theory 
have been observed indirectly b y the continuous ob servation of individual line profiles 



(one-armed prograde precession, 



Telting et al. 



19941 ). and the detection of phase-shifts 
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between different emiss ion line cycles pointing towards a spiral density perturbation 



( IWisniewski et al. 



2007T ). Spectro-interferometric observations of the photocenter shift in 



hydrogen lines of ( Tau appear to directly confirm a on e-armed density mode in the disk 



( iVakili et al. 



1998 



Carciofi et al. 



2009 



StefletaL 



2009|). 



After summarizing the observation and data reduction (Sect. [2]), we will discuss in 
Sect. [3] that our data show that the radial dependence of the hydrogen emission regions can 
be observed directly with KI-SPR, thanks to resolving several emission lines of the Pfund 
series together with Br 7 at unprecedented differential precision in the .fT-band. Concluding 
remarks are given in Sect. HI 



2. Observations 

Details of the observations are given in Table [H Each dataset consists of 155 frames 
of 0.5 sec integration time, taken on the second fringe camera. The frames were selected to 
have an absolute group delay of less then 4.5 /zm to ensure high SNR and li mited impact 



of th e dispersion introduced by the air-filled delay lines (see Appendix C in 



Colavita et al. 



19991 ). The remaining 137 (151) frames of target (calibrator) were corrected for group delay 
(linear) and dispersion (quadratic) wavenumber slopes. The resulting pre-processed frames 
were stacked by calculating the mean in each wavelength bin, and the standard deviation 
of this mean is used throughout the paper as error. This statistical error is adequately 
describing the differential uncertainty of the data. The total flux data were flux normalized 
before the averaging. This is necessary to avoid that the larger absolute flux variation 
due to the spatial filtering of the single-mode fibre fed fringe camera dominates the here 
interesting differential flux variation between adjacent lambda bins. 

The resulting measurables of the SPR observation (flux, V 2 , and dip ) are shown in 
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Fig. [TJ To achieve the rest wavelength calibration shown in the figure, two steps were 
performed. First, the observed wavelength calibration needs to be rectified. The slope 
is measured internally by a Fourier transform based technique. In standard KI observing 
setup, this slope is applied automatically to the data provided through the NExScI archive. 
Furthermore, we used telluric absorption lines in the raw spectra to estimate an offset in 
the wavelength calibration (typically of order of a few tenths of one spectral pixel), and to 
verify the archived rectification. The fact that our data do not show the full K' band is 
due to an alignment problem of the early instrument setup used. For current performance, 
we refer to the KI support webpage Q. In a second step, the corrections for barycentric 



motion of the earth (-13 kins -1 ) and peculiar motion of the star (-6 kms -1 . lEvanslll967l ) 
have been applied. Note for .Br 7 , that the offset between zero dcj) and rest wavelength 
is significantly larger than the accuracy of our wavelength calibration. It is attributed 
to observing a profile which is significantly asymmetric at the level of the used spectral 
resolution. The quality of the resulting wavelength calibration is shown by overplotting the 
Hi-recombination line centers in Fig. [1] (dotted lines). It is reassuring, that (1) the flux 
line profiles are centered, that (2) the various Pfund lines show the same line profile, and 
most importantly, that (3) the dip of both Br 7 and Pfund lines crosses zero close to the 
transition wavelengths, and with the same sign of the slope. (3) is expected for an inclined 
Keplerian disk. 

The target data were calibrated with a subsequently measured unresolved star, 
observed only 8 min later under equal atmospheric conditions. We applied simple 
calibration strategies: the raw target flux was divided by the calibrator, and multiplied 
by a measured spectral template from the IRTF spectral library (A3V, Vacca, personal 
comm.) containing the HI absorption features of the calibrator of similar spectral type. The 



http : //nexsci . caltech.edu/ sof tware/KISupport/v2/v2sensitivity .html 
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spectral template provides a relative flux precision of 1 %, and was folded to R = 1000, 
and divided by a black body continuum fit to the line free regions before application to our 
data. The resultant spectrum (upper panel in Fig. [TJ thus shows in each spectral bin the 
flux with respect to the calibrator. The V 2 was corrected for the system visibility, estimated 
by observing an unresolved calibrator. We did not apply any flux bias correction in our 
data calibration. Although the primary fringe tracker at KI has a flux bias of a few percent 
per magnitude internal tests have demonstrated that the detector used for these science 
observations does not display such a bias. The differential phase is calibrated by subtraction 
of the calibrator dip . To show the impressive differential stability of the data over the used 
band, we overplotted in red linear continua, fitted to the line- free regions (A < 2.15 //m 
and 2.18 < A < 2.315 /im), over flux ratio and V 2 . The slightly red slope of the linear 
flux ratio continuum probably originates in the shell-induced NIR excess of 48 Lib, which 
photosphere is hotter than the calibrator of spectral type A4V. The continuum trend in the 
visibility is due to resolving the innermost circumstellar gas, and is discussed in Sect. 13.11 
The red line in the dip panel marks zero phase, and is not a fit. The resulting pixel-to-pixel 
errors are Afluxratio = 0.003, A V 2 = 0.007, and A dip = 3mrad. The higher relative 
precision of the d(p signal with respect to the precision of the visibility is expected for such 



spectro-interferometric measurements (IMillour et al. 



200J). 



This differential stability at the level of a few 10~ 3 is unprecedented for spectrally 
dispersed OLBI data at the spectral resolution and sensitivity offered by KI-SPR. It is the 
key to derive the scientific results presented in the next section. A more comprehensive 
descriptio n of the SPR mode a nd its sensitivity, performance, and technology will be given 



elsewhere ( Woillez et al. 



20101 ). 



2 http : //nexsci . caltech . edu/ sof tware/KISupport/ dataMemos/f luxbias . pdf 
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3. Results 

Fig. [TJ shows that different components have been detected in the spectra of 48 Lib: 
a clearly resolved continuum emission, and HI recombination line emission of Brj and 
various Pfund lines. The visibility data suggest that the continuum and the line emission 
stem from spatially separated, distinct regions in the circumstellar disk. Also, they derive 
from different emission processes. Therefore, we discuss in the following the two components 
separately. 



3.1. Continuum emission 



The disk continuum emission, emitted by bound-free and free-free thermal emission 
in the disk, is clearly resolved. To facilitate the comparison of our results with other 
published disk diameters, we give the size of a disk-only fit to the continuum visibilities 
(^cont HM = 0-95 ± 0.03 mas) in Table |2j We model the disk continuum by a maximally 



inclined Gau ssian profi 



polarization (IMcDavid 



e, wi th a position angle of 50° which is orthogonal to the axis of 



19991 ) . However the measured continuum emission consists of both 



stellar and disk contributions. Taking into acc ount, that the disk co ntinuum of 48 Lib 



contributes one third of the i^-band emission (iDougherty et al.lll994h . the resulting disk 
size scale (^^a^k = 1-^5 ± 0.05 mas) is seven times larger than the stellar photospheric 
diameter. Source variability may compromise this estimation of fl^^k- But a flux 



calibration of the KI photometry of 48 Lib, based on 



4D 145607 is consistent with the 



2MASS point source flux of our target ( jSkrutskie et al. 



2006) to better than 10%, which is 



the precision level of our absolute flux calibration. iLefevre et al.l (120091 ) report an irregular 
variability in the Hipparcos band at the 5 % level. Therefore, we can expect that variability 
of the star-disk flux ratio has no strong effect on the ^J^^k value reported above. Note 



that although important for the estimation of the disk size, we cannot derive the NIR disk 
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Calibrated SPR data: 48 Lib 
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Fig. 1. — Top: Mean calibrated flux ratio between the 48 Lib and the continuum divided 
calibrator. The red solid line marks a linear continuum fit. The different line profiles of .Br 7 
and P /-emission lines are clearly visible Center: Mean calibrated V 2 of 48 Lib showing that 
both the NIR continuum and the recombination line emission are spatially resolved by the 
interferometer. A linear continuum (red line) was fitted to the line free regions. Bottom: 
The calibrated differential phase data. The red line marks zero phase. All plots show the 
rest wavelength of the target. The vertical dotted lines indicate the rest wavelength of the 
recombination lines. It is apparent that all lines show the same slope at the line center, as 
expected for disk emission. 
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excess from the visibility data due to our small u, f-coverage. Similarly, two-dimensional 
physical or phenomenological disk models of the con tinuum emission are not well constraint 



by a single baseline measurement. 



Gies et al. 



( 120071 ) showed for several Be stars, that, even 



for significantly two-dimensional u, f-coverages, physical models, using radiative transfer 
calculations, and phenomenological models, using Gaussian intensity profiles are difficult to 
distinguish by the visibility data alone. Physical and phe nomenological FW HM estimates 



of the average disk compactness may differ by 50-100 % (IGies et al 



20071 ). This may be 



used as an accuracy estimate for a physical interpretation of the here derived ^^^k 
Since near-infrared disk excess and #^™ k of 48 Lib are within the typical range of Be star 



disks, we can expect that a physical modeling of the disk continuum, as done by 



Gies et al. 



( 120071 ) . would lead to similar densities as described therein. 



The estimated disk continuum-to-stellar radii ratio matches oth er interferometri cally 



observed edge-on systems with comparable properties (e.g. ( Tau in 



Gies et al 



20071 ). The 



typically found moderate extent of the if-ban d emission size (5-10 times the stellar radius) 



is significantly smaller than recently modeled (IStee fc Bittar 



20011 ) . confirming the need for 



interferometric constraints to understand the physical conditions in the stellar outflows of 
Be-stars. 

There is several indications in our data, that the disk continuum emission zone is inside 
the Br 7 and smaller or equal to the Pf— recombination line emission zones. All fitted 
linear characteristics of the violet and red part of the Br 7 (Pf avg ) line emission zones 
(that is 6 I 1 ^ HM and the photocenter shifts PC , for details see the next section; violet 
and red labels negative and positve Doppler shifts, following standard Be-star terminology) 
indicate that the majority of the line emission is created at stellocentric radii larger than 
(comparable to) the radii dominating the continuum emission. Furthermore, the -P/ av g line 
center in total and correlated flux drops below the continuum level (Fig. [3]). This points 
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to a significant absorption of the continuum flux by gas in front of the continuum source, 
that is the continuum emission is embedded and su rrounded by the z one dominating the 
line emission. There 



continuum size) and 



ore, we confirm t he findings of 



Gies et al 



Stefl et al. 



( 120071 . studied Ha versus 



(120091 ) that the observed Bvy (and P/ avg ) hydrogen line 
emission in Be stellar shells appears to occur at larger (or equal) stellocentric radii than the 
disk continuum emission. 



3.2. HI recombination lines 

The emission line widths of Be-stars are dominated by large-scale kinematic Doppler 
broadening (due to the fact that the ionized gas is located in a rotating disk around the 
star) , without significa nt contributions of non-kinematic scattering processes (IHanuschik 



1989 



Gies et al 



20071 ). The disk of 48 Lib is known to be edge on, which simplifies the 



1999 ). 



geometric and kinematic interpretation of the line profiles (IHanuschik fe Vranckenl 
The spectrometer does not fully Nyquist-sample the double peak of the Br 7 and Pfund 
lines, but the central dip in the flux and the shape of the dip signal show that the violet (V) 
and red (R) part of the disk are clearly separated by the spectrometer, and that meaningful 
V/R ratios can be derived. 

P/-24..34 have been detected for the first time in an interferometric spectrum. The 
five strongest individual P /-lines are clear detections (~ 5 — 8 er) above the differential 
precision level of flux ratio, visibility, and d<p . To improve on the SNR for the Pf— lines, 
and to minimize the effects of the limited sampling of the line-profile, we average over the 
five strongest, detected Pf— lines (24.. 28). To do so, the Pf— 25.. 28 lines were scaled to 
match flux levels of Pf— 24. This avoids that the final average profile (-P/avg ) is affected by 
the individual line fluxes. We checked that none of the features of Pf avg discussed below is 
significantly altered or vanishes when averaging over all 11 detected Pfund lines. However, 
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48 Lib: BrGamma-profile 




-1000 -500 500 1000 

vel [km/s] 



Fig. 2. — Zoom into the Br 7 line. The flux ratio (top) is normalized to the linear continuum 
shown in Fig. [TJ The central panel shows the correlated flux, normalized to the linear 
continuum fit. Note that the linear flux ratio continuum fit may overestimate the actually 
observed continuum emission due to absorption by gas in front of the continuum emission, 
which would result in underestimated fluxes at the line center. The measured differential 
phase is plotted in the bottom panel. At the top right of each figure, the integrated violet and 
red line profile is given. The red area marks the region of integration (100-250 km/s), which 
excludes the uncertain central continuum absorption. The wavelength-differential precision 
is indicated by the red errorbar in the left of each panel. 
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Averaged Pfund line-profile 




-1000 -500 500 1000 

vel [km/s] 



Fig. 3. — The three panels show the -P/ av g profile at the flux levels of Pf-(5-24) in the 
same way as in Fig. [2j The smoother line profile with respect to Br 7 derives from the line 
averaging over the linearly interpolated individual Pfund lines. The hatched area indicates 
the beginning of the velocity region where the shape of the mean line profile is significantly 
compromised by the neighboring lines due to the averaging. The profile at large velocities 
(> 330 km s -1 ) are a result of our data processing, and unreal. We show this part here to 
demonstrate that the central, real line profile is more significant than these relatively large 
profile changes in the outer velocity-space compromised from the averaging. The statistical 
error of the mean shapes shown is indicate by the red errorbar in the left of each figure. 
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our choice of averaging over the five strongest lines to create Pf avg leads to the highest 
SNR. In addition, the omitted higher-order Pfund lines are closer in rest wavelength, and 
therefore, overlapping of adjacent line profiles would become apparent and compromise 
the profile of Pf avg . Similarly, at velocities > 330 kms -1 (indicated by the dashed region 
in Fig. [3]), the average profile is affected by adjacent lines. We purposefully show in the 
figure a large velocity space including these artifacts to demonstrate that P f aYg emerges 
even out of these averaging artifacts. Unless indicated otherwise, we use a linear flux ratio 
continuum fit. Using a black body instead of a linear shape for the flux ratio continuum 
leads to a comparably good fit, and slightly increases the off-center line emission of Pf avg • 
We include this uncertainty in the shape of the continuum at the red end of the K'-band in 
the errors given in Table [21 The statistical precision of the averaged profiles are shown as 
red errorbars. 

We compared the individual properties (flux, visibility, and differential phase) of the 
Pf— (24.. 28) included in Pf avg to check if meaningful geometric constraints on the Pfund 
line emission zone can be derived from Pf avg . The decreasing signature in visibility and 
phase with increasing excitation level is consistent with the Pf-n/Pf— 24 line flux ratios. 
This means that a similar geometry is encoded in the measured visibility and phase of 
each line, suggesting a common origin of the lines, and justifying the averaging. The 



individual line ratios used match the respective Pfu nd-line ratio tre nds measured 



Be-star envelopes, such as 7 Cas (IHony et al 



2000) and BK Cam flGranada et al. 



'or other 



2010|). 



Furthermore, these flux ratios vary by no more than up to about 30 % (Pf-28/Pf— 24), 
which indicates together with a similarly moderate variation of the respective Einstein 
coefficients that the lines included in Pf avg are emitted under similar conditions. Therefore, 
they most likely originate from similar regions in the disk. The geometric constraints, 
derived from Pf avg , cannot be more accurate than this level, which matches the numbers, 
given in Table [2j 



- 16 - 



The precision of the V j i?-estimate of the individual Pfund lines is clearly limited 
by our spectral resolution and the sampling of the line profile. Apparent ratios of up to 
20 % above and below unity for the Pfund lines included in P/ avg are regarded as not 
significantly different from unity Therefore, the lines show comparable V/R. Given the 
previous arguments for a common origin of the individual lines in the disk, we can expected 
that averaging does improve the V/P-signal for the Pfund emission zone. 

We concentrate on answering the key question qualitatively: are KI-SPR data suitable 
to measure directly, and efficiently in a single dataset, core parameters like radial extent 
and radius-dependence of the often discussed slowly rotating density inhomogeneity in the 
gas around Be stars? To do so, violet and red velocity bins outside of the line center were 
defined to derive average V-R properties of the lines, while avoiding the zone of highest 
absorption and uncertain continuum (see below). The velocity bins enclose the maximum 
line emission with a width of 150 kms -1 , which is chosen to be comparable to the spectral 
resolution (bins are marked as colored background in Fig. [2^43]). The profile averages for 
each bin are given in the figures. We derived intrinsic line emission sizes and photo-center 
shifts from the data (Table [2]) by assuming the superposition of a linear continuum and 
additional Br 7 dis k emission (for the derivation of the relevant formulae see Appendix C in 



Weigelt et al. 



20071 ) . Note that the intrinsic photocenter shifts are on the order of the fringe 



spacing (X/2B), which requires the use of the exac t complex relation between measured 



and intrinsic phase (Eq. C.3 in 



Weigelt et al. 



20071 ) to derive the correct PC ^ e ft . Table M 



reports the mean linear properties fitted to the linearly interpolated data for each velocity 
bin. 

For our analysis, we used a linear continuum extrapolation of the continuum flux ratio 
into the line center to derive the quantitative constraints in Table El However, the actual 
continuum at the Brj line center may deviate from this linear continuum. The stellar 
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continuum has some photospheric Bvy absorption, and since we observe an edge-on disk 
in emission with, at least partially, optically thick Br 7 emission, additional absorption and 
emission makes it difficult to estimate the exact continuum shape at the line center. This 
presents a difficulty for isolating the correlated flux, visibility and intrinsic dip of the line 
emission at the line center from the measured data, which contain a flux-weighted average 
of both continuum and line emission. An absorbed continuum profile (in contrast to the 
here assumed linear extrapolation of the off-line continuum) at the center of .Br 7 would 
increase the (correlated) line flux, and therefore decrease the derived intrinsic photocenter 
shift. To quantify this, we also applied a B4III .Br 7 template to the data (see also footnote 
( c ) in Tabled]). This would not change the measured V/R asymmetries significantly, but it 
would increase the given 0™ nM by about 40 % and decrease PC j^j? by 20%, which gives 
the order of magnitude of this systematic uncertainty in our quantitative analysis. 



3.3. Radial disk structure 



Br 7 and Pf avg profiles differ. Above all, they show different V/R ratios and central 
absorption levels. The total and correlated .Br 7 line profiles are with V/R ~ 1.8, and 1.6 
respectively, significantly asymmetric. In contrast, ^^ HM a nd PC o f the violet and 



Stefl et al. 



fl2009[ ) resolved the 



red Br 7 emitting regions coincide within the uncertainties. 
Br 7 disk emission of ( Tau with the VLTI with a spectral and spatial resolution similar 
to our data. They found a larger intrinsic line photocenter shift of the brighter wing, 
emitted, at the time of their observation, from the south-eastern, red part of the disk of 
( Tau. A non-LTE disk model can explain such a trend in the dip as well as an extensive 
multi-wa velength and multi- technique dataset with a global oscillation of a spiral density 



20091 ). Our data might show a similar trend of a slightly larger PC 



pattern ( ICarciofi et al. 

1^ of the brighter side of Br 7 . However, this trend is not significant with respect to the 
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Ha (not measured} 



Fig. 4. — Visualization of the measured disk properties (not to scale). The shown one-armed 
over-density pattern would explain the mea s ured V/R and correlated flux profiles of Br 7 



and -P/avg , as predicted by 



Papaloizou et al. 



( 119921 ) , but the here shown pattern is not based 



on a model calculation. The exact shape of the pattern is not constrained by a single KI-SPR 
dataset. The relative location of the optical Ha line is a dded for completeness, matching 



previous, single telescope velocity and V/R measurements ( jWisniewski et al. 



20071). 
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estimated uncertainties. We would encourage a similar, complex modeling effort of the disk 
of 48 Lib, including the new differential constraints from the Pfund emission and based 
on repeated observations to improve on the significance of the d(f> trend of fir 7 . Here we 
concentrate on the overall line emission zone properties. 

The bulk of the fir 7 emission appears to come from similar stello-centric radii, 
and the enhanced emission on the violet side indicates a locally enhanced density. The 
observed asymmetry is inline with the ~9 yr periodicity of cyclic changes of V/R(fia) 



( iHanuschik et al. 



,8 



19951 ). The Ha spectra of 48 Lib, provided by the BeSS database □ also 
report a V/R > 1 in 2008. In contrast, the P/ avg profile appears to be symmetric in all 
properties within the errors. Again, 8™ nM &nd PC f-^ °f P/avg axe comparable for both 
sides, but the smaller PC points to significantly smaller stello-centric radii of the bulk 
of the -P/avg emission, if compared to fir 7 . Both measured V/R(Brj ) and V/R(P/ avg ) 
can be reconciled with a one-armed density perturbation, if a radial dependence of the 
perturbation is allowed, as in the case of a spiral density wave, precessing through the disk. 

There is several indications for that the Pfund emission emerges from inside the bulk 
of the I?r 7 emission. The mean PC (fir 7 ) is with 2.1 mas about twice as large 
as PC jjj^ (-P/avg )• Furthermore, the bins are centered at different stello-centric radii 
and velocities. The peak velocities appear offset (fir 7 : 150 ± 50kms _1 versus -P/avg : 
230 ± 30kms _1 ). These velocity and PC f-^ numbers are consistent with Keplerian gas 
motion, a typical assumption for Be star disk. The Keplerian rotation assumption is further 
supported by the measured double-peaked line profiles (in contrast to highly disordered or 
asymmetric profiles). 

Since Be-star disks are (at least partially) optically thick in the hydrogen emission 



3 http : / /basebe . obspm . f r/basebe/ 
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lines (see the central absorption in Fig. [3] and the Be-star disk study of iLenorzer et al. 
20021 ). it is expected that the stellocentric radius of the respective line anti-correlates 
with the specific line absorption properties. Comparing the respective Einstein absorption 
coefficients, we find -85,24 < -84,7 < -62,3, which suggests that the P/ avg is emitted inside 

of the Br 7 as described by our data. Bvy shoul d be emitted at smaller stellocentric radii 

I 1 F~l 

y con firmed by the data of IWisniewski et al.l (120071 ) (similarly for 



CTau: 


WI , ™ — „~ — 
Carciofi et al. 


, — - 
2009) 


Wisniewski et al. 


2007 


, su 



2009). Also the smaller H a peak velocities (the spectrum in Fig. 1 of 



20071 suggests about 110 km /s) are consistent with larger stellocentric 



radii compared to Br 7 and Pf& vg ■ 

Another indication for the different location of the Pfund line emission are the total 
and correlated flux profiles. Both show a clear absorption below the continuum emission 
level. This, together with the symmetric line profile might indicate that the gas dominating 
the Pfund emission is located in front of the inner disk a nd the star (responsible for 



the continuum) (se e the extensive line profile modeling of 



Hummel fc Hanuschik 



1994 



Telting et al. 



1994J). As discussed, we cannot estimate what fraction of this continuum 
absorption is due to intrinsic photospheric absorption. But a domination of the measured 
continuum absorption by photospheric absorption is unlikely for a Be star. Having the 
majority of the Br'j emission at larger stellocentric radii on the violet side, but the Pfund 



emission at smaller radii in front o 



a radial, one-armed spiral pattern (IWisniewski et al. 



the star suggests that the over-density pat tern follows 



2007 



Carciofi et al. 



20091 ). We depict 



the radial structure of the disk emission as measured by the KI in a schematic way in 
Fig. HI However, our single dataset obviously cannot distingiush a spiral wave from other 
radius- dependent perturbations patterns. 

Our data pr o vide a further quantitative test of Papaloizou's density wave model. 



Papaloizou et al 



(Il992[ ) discuss that the temporal prediction of V/R cycles of about 
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10 yr periods fits the available spectroscopic observations. Interferometric observations 
directly resolve the expected linear scale of the predicted density waves, as opposed to 
indirect methods such as polarimetry (as discussed for instance in IWisniewski et al.ll201( 
The mean photocenter shifts result in stellocentric radii of about 18 (8) R* for Bvy 
(-Pfav g ), which is consistent with both hydro gen emission line modeling of Be stars (e.g. 



Hummel fc Vrancken 



2000; iHonv et al. 



20001 ) and with density wave models predicting 



that the local density perturbation, i nduced by the rotation al flattening of the Be stars 



disappears outside a few stellar radii (jPapaloizou et al. 



1992|). 



4. Conclusions 

We presented the first interferometric detection of several Pfund lines at the red end 
of the K'-band. Thanks to the simultaneous detection of Br 7 in the spectrum of a single 
KI-SPR observation, we can derive stello-centric radius dependent properties the gas disk of 
the classical Be-star 48 Lib at high precision. The geometry and kinematics of the emission 
zones of .Br 7 and of eleven higher Pfund lines are clearly detected in the differential 
visibility and phase signals, and the Pfund emission originates from smaller stello-centric 
radii than .Br 7 . The radial separation of the various emission lines is convincingly 
explained by the physics of optically thick line emission in the Be circumstellar disks. In 
addition, the if'-band disk continuum emission was resolved as well, and shown to be 
emitted from inside the Br 7 (V/R > 1) hydrogen line zone, at radii comparable to the bulk 
of the P/-emission (with V/R ~ 1). Our findings match qualitatively and quantitatively 
theoretical Be-star disk model calculations, predicting a precessing one-armed over-density 
pattern. The spatial constraints from our interferometric data coincide with the common 
notion that Be-star disks have a Keplerian-like radial velocity profile. Observations similar 
to the here presented will shed more light on the enigmas of the creation and properties of 
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circumstellar Be shells at all stages of evolution. In particular time-resolved monitoring of 
Be circumstellar environments with KI-SPR together with radiative transfer modeling of the 
disk emission will deliver strong constraints on the actual radial and azimuthal geometry of 
the disk inhomogeneities likely to be responsible for the HI line emission asymmetries. 

In addition to the presented scientific results on the disk of 48 Lib, we demonstrate 
with the detection of the Pfund lines, that even weak line emission can produce a significant 
signal in the differential phase. KI-SPR data of bright stars can be calibrated to a 
differential precision of Afluxratio = 0.003, A V 2 = 0.007, and Ad(f) = 3mrad. This 
extreme precision in the differential phase equals an on-sky centroid shift of 3 //as (a relative 
precision of 10~ 3 ), which is an improvement by two orders of magnitude over state-of-the-art 
single-telescope spectro-astrometry. This improvement is however limited to the relatively 
small interferometric field of view (~30 mas) and sensitivity of the current KI-SPR mode. 
Dual field phase referencing, the second phase of the KI-ASTRA upgrade, is currently 
under construction, and will allow fainter science targets to be observed. Using a different 
dispersive element, even higher spectral resolution could be achieved in SPR mode. 

We are grateful to W. Vacca and A. Seifahrt for many useful discussions and 
contributions to this work. The excellent support of the KI team at WMKO and NExScI 
helped making these observations a success. The realization of the KI-ASTRA upgrade is 
supported by the NSF MRI grant, AST-0619965. The data presented herein were obtained 
at the W.M. Keck Observatory, which is operated as a scientific partnership among the 
California Institute of Technology, the University of California and the National Aeronautics 
and Space Administration. The Observatory was made possible by the generous financial 
support of the W.M. Keck Foundation. The authors wish to recognize and acknowledge 
the very significant cultural role and reverence that the summit of Mauna Kea has always 
had within the indigenous Hawaiian community. We are most fortunate to have the 



-23 - 



opportunity to conduct observations from this mountain. The Keck Interferometer is 
funded by the National Aeronautics and Space Administration as part of its Exoplanet 
Exploration program. This research has made use of the SIMBAD database, operated at 
CDS, Strasbourg, France. This research has made use of NASA's Astrophysics Data System 
Bibliographic Services. This work has made use of the BeSS database, operated at GEPI, 
Observatoire de Meudon, France. 

Facilities: Keck:I (), Keck:II () 
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Table 1: Observing log from April 25, 2008 (UT). 



Name 


v/h/km 


UT BL pro j 


Stellar diam 








[mas] w 


48 Lib (Tar, B4IIIW) 


4.9/4.8/4.6 


13:20:49 66.6 m, 39° EoN 


0.23 ± 0.1 


HD 145607 (Cal, A4V) 


5.4/5.3/5.1 


13:28:17 


0.25 ± 0.04 



"From Simbad 

b Bolometric photosphere fit using NExScI's getCal tool. Diameters < 0.5 mas are unresolved by the KI. 
c From 



by 



Floauet et al 



Slettcbakjj 19821). However the spectral type is uncertain by a few sub-classes as discussed for instance 



1996). 
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Table 2: Measured and derived properties of the circumstellar shell of 48 Lib. 



^cont (^aUSS- 


■disk 


only): 


ziFWHM 
^cont 


n n a _i_n no rt ^, 
U.y4±U.Uo mas 




^cont (Gauss- 


■disk 


+ star): 


/DFWHM 
^cont^isk 


1.65±0.05 mas 




Pr 7 






^(-250,-100 kms- 1 ) 


P(100,250 kms- 1 ) 


V 1 R 


Fii n e (norm.) 






0.29 ±0.01 


0.16 ±0.01 


1.80 ±0.2 


(norm.) 






0.23 ±0.01 


0.14 ±0.01 


1.66 ±0.2 


line (mas) 






1.7 ±0.2 


1.4 ±0.2 


1.2 ±0.2 


PC ^ (mas; 


) 




2.2 ±0.2 


1.9 ±0.2 


1.2 ±0.2 


P/av g (5 : 24. 


.28) 




^(-325 , -175 kms- 1 ) 


P(175,325 kms- 1 ) 


V/R 


F line (norm.) 






0.029 ±0.01 


0.028 ±0.01 


1.0 ±0.1 


F£°™ (norm.) 






0.024 ±0.01 


0.024 ±0.01 


1.0 ±0.1 


C™ (mas) 






1.9 ±0.3 


1.5 ±0.2 


1.2 ±0.3 


PC (mas; 


) 




0.9 ±0.3 


0.9 ±0.3 


1.1 ±0.3 



a All uncertainties given in the table are derived from the standard deviations of the mean total flux and 
correlated flux spectra shown in the figures. These uncertainties reflect the differential, relative precision of 
the values, and is adequately used throughout the article for intcrcomparison of the different components. 
However, the absolute accuracy of the individual estimates is 10% for continuum and Br 7 values, and 25% 
for the P/avg values. 
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